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ABSTRACT 


An  approximate  solution  of  the  diffusion  equation  describ¬ 
ing  electron  diffusion  in  a  multimode  rectangular  cavity  is  described. 
The  results  indicate  that  the  electron  concentration  builds  up  in 
the  spaces  between  the  maxima  of  ionization.  This  means  that  the 
breakdown  electric  field  is  characteristic  of  a  diffsuion  length 
approximating  that  of  the  lowest  diffusion  mode. 
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ELECTRON  DIFFUSION  IN  A  MULTIMODE  CAVITY 


There  have  been  several  studies  of  electron  diffusion  in 
cavities  in  which  the  non-uniform  nature  of  the  microwave  electric 

1  2 

field  has  complicated  our  understanding  of  the  breakdown  phenomena.  9 

3 

In  a  previous  report  to  the  Rome  Air  Development  Center  an  analysis 
of  the  problem  was  given  using  a  variational  treatment.  The  results 
were  presented  in  the  form  of  the  ratio  of  an  effective  characteristic 
diffusion  length,  A  ,  to  the  geometrical  diffusion  length,  A  ,  for 
different  values  of  E/p,  where  E  is  the  electric  field  at  its  maximum 
and  p  is  the  pressure.  We  have  measured  the  breakdown  field  in  a 
cavity  of  the  type  described  in  this  report  and  find  that  experiment 
and  theory  do  not  agree.  In  this  report  we  present  a  simpler  theory 
of  the  breakdown  phenomena  and  the  results  of  the  experiments,  as 
well  as  the  conclusions  which  we  may  draw  from  them. 

We  define  the  particle  current  density,  f  ,  by  the  equation 

r  =  -  \7(Dn),  (1) 

where  D  is  the  diffusion  coefficient  and  n  is  the  electron  concentra¬ 
tion.  The  continuity  equation  is 

+  v  -  r  +  p  =  o,  (2) 

where  P  is  the  production  or  loss  rate  at  sources  or  sinks  within  the 
volume  considered.  We  may  consider  breakdown  to  take  place  when  the 


2 

production  rate  slightly  exceeds  the  loss  rate  and  <9n/*Pt  may  be 
set  equal  to  zero. 

Let  Dn  -  Y  and  combine  Equations  (1)  and  (2).  The  result 
is 

V2  V  +  p  =  o.  (3) 

The  production  rate  within  the  cavity  is  determined  by 

ionization  and  is,  in  fact,  nv^,  where  is  the  number  of  ionizations 

2 

per  second  per  electron.  nv^  is  then  equal  to  £ E  Y ,  where  £ 

2 

is  the  high  frequency  ionization  coefficient,  defined  as  v^/DE 
(E  is  the  rms  value  of  the  electric  field). 

Equation  (3)  then  becomes 

V2  If  +  /e2$/  =  0.  (4) 

The  electron  concentration  vanishes  on  the  walls  of  the 
cavity  and  this  provides  us  with  the  necessary  boundary  conditions 
on  Equation  (4) . 

The  cavity  in  which  we  are  particularly  interested  is  a 
rectangular  parallelepiped  resonant  in  the  Th^Q  mode.  The  electric 
field  is  then  given  by  the  following  expressions,  where  x,  y,  and  z 
are  the  coordinates  indicated  in  Figure  1. 
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FIGURE  1 

SCHEMATIC  REPRESENTATION  OF  TM23Q  CAVITY 

E  =  0  E  =  0 

x  y 

(5) 

E  *  E  sin  — —■  sin 
z  o  A  B 

We  use  the  approximation  used  by  Herlin  and  Brown  for  this 
ionization  coefficient,  namely, 

/  -  /0  (E/E0)P'2,  (6) 

where  Eq  is  the  field  at  the  center  of  the  cavity.  Therefore,  Equa¬ 
tion  (4)  becomes 

V2  Y  +  (sin  Zf-  sin  V  »  o.  (7) 

Typical  values  of  p  range  from  about  8  to  16.  A  plot  of 
si  A  for  two  or  three  values  of  p  is  shown  in  Figure  2.  This  gives 
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an  indication  of  what  a  cross  section  of  an  ionization  surface  over 
the  area  of  the  cavity  would  look  like.  Inspection  shows  that  the 
ionization  is  less  than  one-third  of  the  maximum  over  a  distance  of 
two-thirds  of  the  x  axis.  As  a  first  approximation  we  shall  use  a 
step  function  to  represent  the  ionization,  and  this  is  represented 
analytically  by  the  following  function: 

2 

jfE  =  P,  a  constant  .7a  <  x  <  1.3a 

2.7a  <  x  <  3.3a 

.7b  <  y  <  1.3b  (8) 

2.7b  <  y  <  3.3b 

4.7b  <  y  <  5.3b 

where  A  -  4a  and  B  =  6b  and  is  zero  elsewhere. 

The  dotted  lines  in  Figure  2  indicate  how  the  approximation 
compares  with  the  actual  function. 


21  -x 

1.3  a 
0,1  CL 


* -  - - > 


FIGURE  3 

CROSS-SECTIONAL  REPRESENTATION  CF  APPROXIMATION  TO  IONIZATION  RATES 
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We  can  also  represent  this  ionization  function  as  viewed 
from  the  direction  of  the  z  axis  in  Figure  3.  In  this  figure  the 
cross-hatched  areas  represent  those  areas  in  which  there  is  ioniza¬ 
tion. 

Equation  (7)  then  becomes 

V1  1/  +  P  v  =  o,  (9) 


where  P  *  0  for  most  of  the  plane  of  Figure  3  and  is  a  constant  in 
the  cross-hatched  areas. 

Because  there  is  no  variation  of  field  in  the  z  direction 
we  may  assume  that  the  z  variation  of  Y  is  sinusoidal. 

We  let  Y  =  g(x,y)  sin  ,  where  g  is  a  function  of  x  and 
y  but  not  of  z. 


1 

g 


2 

-y  +  P  =  0. 

C 


(10) 


The  solution  of  this  equation  with  the  type  of  ionization 
depicted  in  Figure  3  is  extremely  difficult  because  the  equation  may 
not  be  further  separable  and  because  of  the  nature  of  the  boundary 
conditions. 

However,  we  can  find  out  a  good  deal  about  the  nature  of 
the  physical  processes  going  on  by  reducing  the  problem  to  one  in 
which  we  assume  that  the  ionization  is  a  constant  equal  to  P  in  the 
areas  .7a  <  x  <  1.3a  and  2.7a  <  x  <  3.3a,  and  zero  elsewhere.  This 
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is,  of  course,  not  very  close  to  what  happens;  but  insofar  as  we  are 
dealing  with  variations  in  the  electron  concentration  in  the  x  direc¬ 
tion,  it  is  a  good  approximation,  at  least  for  some  values  of  y.  We 
can  then  find  the  manner  of  the  variation  in  the  x  direction  and  infer 
that  the  y  variation  is  similar  functionally. 

We  get  g(x,y)  =  X(x)  Y(y),  and  because  there  is  no  variation 
in  the  ionization  in  the  Y  direction  we  set  Y  =  sin  .  Equation  (10) 
then  becomes 

~M"£2'£2  +  P  =  0-  (H) 

X  dx  B  C 

2  2 

When  P  is  greater  than  -y  +  ,  X  is  sinusoidal;  and  when 

BZ  C 

P  is  zero,  X  is  exponential.  Therefore,  the  solutions  are  of  the 
fora 

X  =  cos  (^x  -  <f), 
where  there  is  ionization,  and 

X  =  sinh  ox  4*  k  cosh  ox  elsewhere. 


These  solutions  lead  to  two  conditions: 


and 


2 


a 


2 

n 


2 

a  . 


(12) 


(13) 


We  now  write  down  the  details  of  the  x  variation: 
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x  <  .7a 

X  « 

d^(sinh  ox  +  k^  cosh  ox) 

(14) 

.7a  <  x  <  1.3a 

X  * 

COS  [\±  (x 

-  a)  -  <f] 

(15) 

1.3a  <  x  <  2.7a 

X  = 

^2  j^os k 

[a(x  -  2a)  ] 

+  l^sinh  a(x  - 

(16) 

2.7a  <  x  <  3.3a 

X  = 

COS  |ji(x 

-  3a)  +  cTj 

(17) 

3.3a  <  x  <  4.0a 

X  * 

d^  £sinh  a(x  -  4a)  + 

cosh  a(x  -  4a)J 

(18) 

This  function  must  be  continuous  and  have  a  continuous  first 
derivative  and  it  must  vanish  at  the  boundaries.  Therefore, 

a  =  0.  Because  of  the  symmetry  the  slope  at  x  =  2a,  the  center 
of  the  x  variation,  is  zero  and,  therefore,  k^  is  zero.  These  boundary 
conditions  lead  to  the  following  equations: 


d^  sinh  .7aa  = 

cos(-.3afi  -  oT) 

(19) 

d^a  cosh  .7aa  = 

-  [i  sin  (- . 3a^  - 

(20) 

d^  cosh  (-.7aa) 

»  cos  (,3ap.  -  cf) 

(21) 

d^a  sinh  (-.7aa) 

=  -  p.  sin  (.3ap  -  S ) 

(22) 

The  dimensions  of  the  cavity  studied  are  A  »  4.91  cm, 

B  «  6.36  cm,  and  C  =»  6.98  cm.  The  value  of  the  characteristic  diffusion 
length  for  the  lowest  mode  would  then  be  A  =  1.08  cm.  From  Equation 
(12)  we  find  that  a  s  .6683/cm.  Because  we  are  dealing  with  a  non- 
uniform  electric  field  we  are  unable  to  specify  P  initially.  P  is 
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determined  by  the  value  of  the  high  frequency  ionization  coefficient, 

2 

£ o,  and  the  value  of  the  breakdown  field,  ,  is  equal  to  v/D, 

2 

which  is  equal  to  1/A  in  the  uniform  field  case  and  may  be  set 
2 

equal  to  1/ A  *  where  A  is  some  effective  characteristic  diffusion 
e  e 

length  in  this  case.  As  a  first  approximation  to  A^  we  use  the 
uniform  field  lowest  mode  A  and  then  solve  the  problem  by  successive 
approximations,  finding  values  of  | a,  d^,  d^,  and  cf  which  are  consis¬ 
tent  with  Equations  (19)  -  (22).  The  values  so  obtained  are  as 
follows: 

dl  =  1.20 

d2  »  .830 

\i  =  1.37 /cm 

-  .253 

The  electron  concentration  is  plotted  as  a  function  of  x  for 
these  values  in  Figure  4. 

On  Figure  4  we  have  also  plotted  a  sinusoidal  curve  with 
the  same  peak  concentration. 

We  may  now  reverse  the  procedure  and  calculate  an  electron 
concentration  profile  along  the  y  direction,  and  again  calculate  the 
various  constants  by  successive  approximations.  The  results  which  are 
not  reproduced  here  are  similar  to  those  in  the  x  direction,  and  the 
electron  concentration  profile  has  three  peaks.  The  concentration 
gradient  at  the  edge,  which  is  the  only  number  we  shall  use,  is 
.782/cm.  This  value  and  the  corresponding  value  for  the  x  direction 


10 


FIGURE  4 

ELECTRON  CONCENTRATION  FROM  EQUATIONS  (14)  THROUGH  (18) 
COMPARED  WITH  sin  x 
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are  not  sensitive  to  the  actual  size  of  the  volume  in  which  ioniza¬ 
tion  occurs.  Variations  of  factors  of  four  in  the  size  of  the  volume 
make  very  small  changes  in  the  result. 

We  now  use  these  results  to  calculate  an  effective  diffu¬ 
sion  length.  The  rate  of  electron  loss  at  any  surface  is  determined 
by  the  concentration  gradient,  and  one  way  of  calculating  the  effective 
characteristic  diffusion  loss  would  be  to  calculate  the  total  electron 
losses  at  all  the  surfaces  and  average  the  ionization  rate  over  the 
whole  cavity  volume  and  then  set  the  loss  rate  equal  to  the  ioniza¬ 
tion  rate.  However,  we  may  approximate  this  result  closely  by  a  much 
simpler  procedure. 

The  ionization  rate  is  proportional  to  the  electric  field 
raised  to  a  very  high  power  and,  therefore,  since  the  field  need 
only  change  very  little  to  make  a  large  change  in  ionization  rate, 
we  will  not  need  to  average  over  the  cavity  volume.  The  characteristic 
diffusion  length  is  ordinarily  defined  by  the  equation 

l/A  =  Yi  +  Y2  +  Y3  > 

where  y V2*  and  Y3  are  the  eigen  values  of  the  function  describing 
the  electron  concentration.  In  the  lowest  diffusion  mode  case  they 
also  represent  the  slope  of  the  electron  concentration  at  the  central 
point  of  each  surface;  i.e.,  the  concentration  is  n  *  nQ  sin  y^x 
and  *  \y ^  cos  yx~\  We  can,  therefore,  make  a  good 

approximation  to  the  value  of  the  effective  diffusion  length  by 
determining  the  three  y'8  by  finding  -77^  at  the  center  of  each  surface. 

Cr  X 
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Inspection  of  Equation  (14)  gives  us  y^  *  d^a  31  1.00,  and 
the  corresponding  result  for  the  y  direction  is  y^  =  p  =  .782.  There¬ 
fore,  we  find 

A  =  .743  cm. 

e 

Using  this  value  of  Ag,  we  calculate  the  high  frequency 
ionization  coefficient,  and  plot  it  as  a  function  of  E/p  for 
various  values  of  p\,  along  with  the  values  obtained  for  air  in  flat 
cylindrical  cavities  of  dimensions  such  that  the  electric  field  is 
uniform.  This  is  shown  in  Figure  5,  and  the  original  data  from  which 
£  is  calculated  are  shown  in  Figure  6.  Also  shown  in  Figure  6  are 
breakdown  fields  in  the  multimode  cavity  for  nitrogen  and  hydrogen. 

The  nature  of  the  resulting  curves  shows  that  the  results  we  obtained 
are  reasonable,  although  at  very  high  E/p  the  value  of  Afi  appears 
definitely  too  low.  High  E/p  corresponds  to  low  pressures  and  at  the 
low  pressures  it  is  likely  that  the  electron  concentration  profile  is 
closer  to  that  of  the  lowest  mode,  and  this  would  lead  to  a  smaller 
value  of  /  .  The  results  obtained  for  the  variational  calculations, 
on  the  other  hand,  do  not  all  agree  with  the  results  from  the  uniform 
field  experiments. 

The  variational  calculations  were  made  with  the  assumption 
that  the  electron  concentration  was  very  low  in  the  region  between  the 
ionization  peaks,  whereas  the  present  calculation  assumes,  as  we  see 
from  Figure  4,  that  the  regions  between  high  ionization  points  fill 
up  with  electrons  because  there  is  no  place  to  which  they  can  diffuse. 


C  (VOLTS) 
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FIGURE  5 

HIGH  FREQUENCY  IONIZATION  COEFFICIENT  /  AS  A  FUNCTION  OF  E/p 
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IN  RECTANGULAR  TM,,_  CAVITY  AT  9.4  KMC 
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